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ABSTE5CT 

A  coupled  oceanic-atmospheric  boundary  layer  model  which 
provides  single-station  prediction  capability  is  evaluated 
relative  to  wind  stress  and  oceanic  mixed  layer  depth  obser- 
vations. The  model  is  initialized  and  verified  using  data 
obtained  during  the  1983  Mixed  Layer  Dynamics  Experiment 
(MILDEX).  Model  computation  of  friction  velocity  (u*)  and 
oceanic  mixed  layer  depth  are  compared  with  observations 
from  both  atmospheric  frontal  and  non-frontal  synoptic  situ- 
ations. Favorable  results  of  model  u*  predictions  are 
achieved  although  in  some  cases  they  are  slightly  higher 
than  observed  turbulence  derived  values  of  u*.  Mixed  layer 
predictions  are  very  close  to  the  observed  except  in  frontal 
regions.  Proper  ABI/OBL  model  initialization  of  boundary 
layer  values  is  critical  in  order  to  achieve  favorable 
results. 
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I-  INTEODECTION 

Naval  oceanographers  and  meterologists  obtain  a  great 
variety  cf  informaticn  to  help  them  in  making  atmospheric 
and  oceanic  forecasts.  Extensive  use  of  large  shore  tased 
computers,  remote  sensing  eguifment  and  world-wide  communi- 
cations networks  make  the  forecaster  rely  on  remotely  tased 
technology  rather  than  local  indicators  and  have,  to  a  great 
degree,  removed  the  "art"  from  single-station  forecasting. 

However,  large-scale  conflict  or  operational  considera- 
tions could  deny  all  forecasting  inputs  available  to  the 
forecaster  except  local  indications.  Reliance  on  single- 
station  forecasting  for  naval  operations  during  times  of 
crisis  or  war  was  first  addressed  by  Oliver  and  Oliver 
(1945): 


During  the  last  few  years,  wartime  conditions  have  made 
it  necessary  for  isolated  combat  units  to  issue  fore- 
casts in  regions  where  no  network  of  meteorological 
stations  could  be  available.  Frequently  the  data  from 
several  stations  or  from  reconnaissance  planes  are 
available,  but  in  seme  regions  the  forecaster  must  rely 
only  on  surface  and  upper-air  observations  made  at  his 
own  station.  This  is  particularly  true  in  the  case  of 
ships  at  sea.  Hence,  it  is  important  to  develop  profi- 
ciency in  extracting  information  from  limited  aerolo- 
gical  data. 


While  the  above  quote  pertains  to  atmospheric  forecasts,  the 
same  could  be  true  for  ocean  predictions. 

In  modern  warfare,  isolated  ships  or  battle  groups  far 
from  the  main  force,  operating  in  an  emission  control 
(EMCON)  posture,  must  be  able  to  predict  their  operating 
environment  without  shore  support.  Such  predictions  are 
important  for  the  effective  utilization  of  weapon  and  sensor 
systems  and  to  deny  the  effective  use  of  the  environment  to 
the  opposing  force-    Affected  systems  include:    navigation 
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radars,      communications      equipment,       over-the-horizon       (OTH) 
radars,    air   and   surface   search    ladars,      missile   and    gun   fire 
control   radars,      laser  and      infra-red    (IR)       guided    missiles, 
hull    mounted  sonars,    towed   arrays,    and    variable   depth    sonars 
(YDS)  . 
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Fig.  1-1    Components  of  the  Coupled  Atmospheric 
and  Oceanic  Boundary  Layer  Model. 


A  coupled  atmospheric  and  oceanic  model  (O'Loughlin, 
1982)  has  been  developed  to  predict  ocean  and  atmospheric 
parameters  from  local  indications.  The  components  of  this 
coupled  model  are  depicted  in  Fig.  1.1.   This  model  is  being 
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tested  and  evaluated  with  shipboard  applications  in  mind. 
Accurate  predictions  of  ocean  mixed  layer  depths  (MLD) , 
thermocline  jumps  and  gradients  are  essential  in  optimizing 
anti-submarine  warfare  (ASff)  systems  and  employing  acoustic 
count ermeasures.  These  parameters  are  dependent  on  the  heat 
fluxes  and  wind  stress  at  the  ccean  surface.  This  thesis 
will  examine  how  the  single-station  coupled  model 
predictions  of  the  surface  fluxes  and  mixed  layer  depth 
compare  with  observations. 
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II.  MODEL  BACKGBODND  ANE  ANAUSIS  PBOCEDOBES 

A.   BOONDABT  LAIEB  DESCBIPTION 

Thermodynamically  and  dynaaically  forced  marine  atmos- 
pheric (MABL)  and  oceanic  (OBI)  boundary  layers  must  be 
considered  in  coupled  model  description.  The  air-sea  inter- 
face has  adjacent  oceanic  and  atmospheric  turbulent  mixed 
layers  which  effectively  insulate  the  guasi-geostrophic 
oceanic  and  atmospheric  flows  atove  the  top  of  the  MABL  and 
below  the  bottom  of  the  OBL.  The  primary  source  of  the 
turbulence  is  the  velocity  (shear)  and  buoyancy  (density) 
gradients  created  by  the  exchange  of  mass,  energy  and 
momentum  across  the  air-sea  interface.  Such  turbulence 
generates  turbulent  kinetic  energy  (TKE)  and  vigorous 
vertical  mixing  (Davidson  and  Garwood,  1984). 

Under  undisturbed  conditions,  the  MABL  consists  of  a 
moists  well-mixed  layer  extending  from  the  surface  to  the 
capping  inversion.  At  the  inversion,  jumps  in  temperature 
and  humidity  occur  with  a  rapid  increase  in  temperature  and 
a  corresponding  decrease  in  humidity.  Within  the  MABL, 
equivalent  potential  temperature  and  specific  humidity  are 
conservative  quantities. 

A  consequence  of  the  large  velocity  fluctuations  and 
mixing  in  the  MABL  is  that  energetic  eddies  (extending  from 
the  surface  to  the  inversion)  entrain  warm,  dry  air  and 
bring  momentum  into  the  mixed  layer.  This  results  in  the 
upward  growth  of  the  layer.  Subsidence  arises  from  large- 
scale  atmospheric  forcing  and  tends  to  limit  the  depth  of 
the  mixed  layer.  Therefore,  the  change  in  the  inversion 
height  with  respect  to  time  is  a  function  of  the  entrainment 
rate  and   subsidence.    If   entrainment  causes   the  MABI   to 
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extend  above  the  lifting  condensation  level  (LCL) ,  clouds  or 
fog  will  form  within  the  layer  and  greatly  affect  the  OBL- 

In  the  ocean,  the  OBL  is  defined  as  extending  from  the 
surface  to  the  top  of  the  seasonal  thermocline.  Within  the 
mixed  layer,  nearly  homogeneous  profiles  of  temperature  and 
density  are  evident.  While  density  is  a  function  of  both 
salinity  and  temperature,  over  short  time  scales,  tempera- 
ture has  a  far  greater  effect  (Miller,  1976),  At  the  bottom 
of  the  mixed  layer,  large  changes  in  temperature  and 
salinity  exist  with  increasing  depth.  This  is  the  top  of 
the  thermocline  region  and  serves  as  a  transition  from  the 
mixed  layer  to  the  usually  dynamically  stable  interior 
ocean.  Idealized  model  profiles  for  both  the  MABL  and  OBL 
are  depicted  in  Fig.  2.1. 

Radiation  is  another  source  of  energy  in  both  mixed 
layers.  While  most  solar  energy  will  penetrate  the  MASL, 
this  is  not  true  of  the  OBL.  Long  wave  radiation  is 
absorbed  and  emitted  by  the  first  few  millimeters  of  sea 
water.  Therefore  downward  turbulent  heat  flux  is  as  impor- 
tant as  the  upward  flux.  Surface  heat  flux  in  the  OBL  is  a 
function  of  sensible  and  latent  heat  as  well  as  the  long  and 
short  wave  solar  radiation-  This  radiation  is  highly  depen- 
dent on  cloud  formation  in  the  ABL  so  proper  ABL  modelling 
is  essential  for  realistic  ocean  MLD  predictions. 

While  buoyancy  forces  create  turbulence  in  both  boundary 
layers,  relatively  more  turbulence  in  the  OBL  is  produced  by 
wind  stress  (DavidscD  and  Garwood,  1984) .  The  MLD  will 
deepen  if  the  wind  generated  downward  flux  of  momentum  and 
forced  ccnvective  mixing  {of  surface  water  that  has  been 
cooled  by  radiative  heat  loss  or  evaporation)  provide  suffi- 
cient TKE  to  erode  the  thermocline  and  entrain  cooler  water 
downward-  If,  on  the  other  hard,  the  near  surface  absorp- 
tion of  solar  radiation  produces  a  downward  buoyancy  flux 
which  dampens  turbulence  and  provides   for  net  warming,   the 
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Fig-  2.1    Idealized  Atnospiieric  and  Oceanic 
Boundary  Layer  Temperature  Profile. 


MLD  will  shallow.  Fig.  2.2  shews  the  importance  of  mechan- 
ical mixing  to  the  CBL  and  its  interaction  with  buoyancy 
fluxes-  A  primary  objective  of  this  thesis  will  be  to 
determine  the  effects  of  wind  stress  on  the  MLD- 

Several  cause  and  effect  relationships  are  evident. 
First,  clouds  can  be  caused  by  changes  in  the  ocean  surface 
temperatures  which,  in  turn,  affect  the  radiation  budget  of 
the  OBL  which  results  in  surface  temperature  changes- 
Secondly,  the  stability  of  the  iHABL  is  influenced  by  air- 
ocean  temperature  differences.  Again,  cloud  formation  plays 
an  important  role.  Such  relationships  make  a  coupled 
prediction  approach  necessary-  This  coupled  model  uses 
local  oceanic  descriptions   from  a  bulk  oceanic   mixed  layer 
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Fig.  2.2    Bechanical  Energy  Budget  for  the 
Ocean  Hiied  Layer. 


model  and  the  atmospheric  descriptions  come  from  an  atmos- 
pheric boundary  layer  model.  Such  a  coupled  model  has 
tactically  significant  implications  to  forecast  parameters 
that  are  important  to  weapons  systems.  Fig.  2.3  shows  the 
interrelationships  between  the  coupled  model  and  tactical 
products  which  can  use  model  outputs  (on  the  extreme  right 
hand  side) . 
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B.   flAEIHE  ATHOSPHEEIC  BOnNDABY  LATER  (MABL)  MODEL 

Existing  integrated  MABL  irodels  are  zero-order,  two 
layer,  mixed  layer  models  consisting  of  a  well  mixed  turbu- 
lent boundary  layer  underneath  a  relatively  non- turbulent 
free  atmosphere.  A  transitior  zone  (inversion)  of  zero 
thickness  (hence  the  term  zerc-order)  separates  the  two 
layers.  Because  it  has  no  thickness,  a  jump  occurs  at  the 
inversion  of  the  conservative  parameters.  The  model  was 
described  by  Stage  and  Businger  (1981),  who  formulated 
entrainment  energetics,  and  modified  by  Davidson,  et  al. 
(1984),  who  described  bulk  aerodynamic  formulae  for  the 
surface  layer. 

The  following  inputs  are  required  by  the  model; 

(1)  an  initial  atmospheric  sounding; 

(2)  the  mean  winds  at  a  level  within  the  surface  layer; 

(3)  the  surface  temperature;  and 

(4)  the  subsidence  rate. 

As  formulated,  ten  winds  can  be  input  over  the  forecast 
period.  The  sea-surface  temperature  (SST)  remains  unchanged 
and  the  surface  current  is  assumed  to  be  zero.  SST  and  the 
surface  current  change  in  the  coupled  version  due  to  the 
interaction  with  the  OBL  model.  Inversion  height,  mixed 
layer  values  of  temperature,  hunidity  and  wind  are  predicted 
at  30-minute  intervals  in  both  the  coupled  and  uncoupled 
versions-  The  prediction  stefs  are  shown  in  Fig.  2.4. 
Procedures  are  the  same  for  clear  and  cloudy  cases  except 
for  entrainment  computation  and  estimation  of  cloud  top 
cooling.  Formation  of  clouds  cr  fog,  cloud  top  cooling  and 
associated  entrainment  are  important  in  the  physical 
processes   in   the    MABL- 
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Fig.    2.4        Scheaatic   of   Irput,    Prescription   and 
Computing  Steps    in    MaBL    Prediction- 


Tie  MABL  model  is  very  sensitive  to  the  prescriied  large 
scale  subsidence  which  moderates  the  depth  of  the  moist 
marine  layer  and  hence,  cloud  formation.  Methods  which  can 
be  used  to  compute  the  subsidence  from  single  station  obser- 
vations   are:         the    kinematic      method;      the    adiabatic    method; 
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and   the  integratioc   of  the   mcisture   budget  equation   (Q) 
method.   Gleason  (19  82)  found  that  the  Q-method  has  the  most 
merit  and  was   used  in  this  thesis.     However ^   this  method 
assumes  that  no  advection  will  take  place. 

Integrated   rate  equations  by   Tennekes  and   Dreidonks 

(1981)   are  used   to  predict  the  time  rate  of  change  of  the 

conservative  quantities   and  their   respective  jumps   at  the 

inversion.     The   vertically   integrated   equations   fcr   a 

conserved  property,  X,  are: 


h(DX/Dt)  =  (w'x»)  -(w'x'),  +  scurce  (2.1) 

o         n 


h(DAX/Dt)  =  h^x(^h/^t)  -  (w'x»)  +  (w»x»).  -  source   (2.2) 

o  h 

where : 

source  =  -  (F  .  -F  )/pCa      for  X  =  temperature 

source  =  0  for  X  =  humidity 

and  Gamma  (r)  is  the  conservative  property  lapse  rate  above 
the  inversion  and  Fp  is  the  net  radiative  flux.  The 
subscripts  "h"  and  "o"  refer  to  inversion  height  and  surface 
values  respectively. 

The  entrainment  velocity  parameterization  is  based  on 
the  formulation  by  Stage  and  Businger  (198  1)  and  is  used  to 
close  the  system  of  equations  and  determine  the  time  evolu- 
tion of  the  inversion  height.  Ihe  closure  assumpton  is  that 
the  dissipation  rate  of  the  TKE  is  a  fixed  fraction  (1-A)  of 
the  production  rate.  "A"  is  the  entrainment  coefficient  and 
is  taken  to  be  0.2. 

Bulk  aerodynamic  formulas  described  by  Davidson,  et  al. 
(1984),  are  used  to  determine  the  surface  fluxes  of 
momentum,  sensible  heat  and  latent  heat.   They  are: 

u*  =  C^^/^n^Q         (friction  velocity)  (2.3) 
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T*    =   C    1/2  (9   -    9)  (temperature    parameter)  (2.4) 

q*   =   C    1/2  (g  -    q)  (moisture   parameter)  (2.5) 

These   fluxes  are    given   by: 


u»w»    =   u*2  (momentum)  (2.6) 


T'w'    =    u*T*  (sensible    heat)  (2.7) 


q'w«    =   u*q*  (latent   heat)  (2.8) 

where  C      and  C^    are    tlie   stability   dependent    wind    and   temper- 

d  d 

ature  drag  coefficients,  9  is  the  potential  temperature  and 
g  is  the  specific  humidity.  Ihe  subscript  "o"  denotes  the 
surface  value. 

Because  radiation  is  a  primary  factor  in  the  OBI  devel- 
opment, considerable  effort  has  gone  into  calculating  the 
radiation  budget  in  the  MABL.  Oncertainty  in  background 
aerosols,  distribution  and  concentration  of  atmospheric 
absorbing  gases,  and  cloud  droplet  size  spectra  are  sources 
of  error  in  the  MABL  radiation  calculations.  The  short  and 
long  wave  fluxes  are  computed  separately.  Short  wave  radia- 
tive flux  is  calculated  using  the  delta-Eddington  method 
(Joseph,  et  al.,1976)  which  accounts  for  heating  of  the 
mixed  layer  by  solar  radiation.  Incident  flux  at  the  top  of 
the  mixed  layer  is  obtained  fron  the  flux  at  the  top  of  the 
atmosphere  and  the  average  transmittance  in  each  of  15  bands 
covering  the  spectrum  from  0 . 2  um  to  1.7um.  Short  wave 
extinction  is  due  to  both  scattering  and  absorption.  Water 
vapor  in  the  atmosphere  is  the  primary  absorbing  constit- 
uent- Scattered  radiation  forns  a  second  short  wave  radia- 
tive component  (diffuse  solar  radiation)  which,  in  the  mixed 
layer,   is  due  to  atmospheric   particles  (cloud  droplets  and 
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aerosols).  Total  short  wave  radiative  flux  at  the  surface 
is  determined  by  computing  both  the  direct  and  diffuse  radi- 
ative components.  Fairall^  et  al.  (1981)  gives  a  comprehen- 
sive review  of  this  method.  Solar  zenith  angle  is  necessary 
in  order  to  compute  short  wave  radiative  flux  and  it  is 
determined  from  the  latitude,  Julian  day  and  time  of  day  of 
the  initial  input.  The  fraction  of  short  wave  radiation 
from  the  sea  surface,  Ag,  which  is  an  important  short  wave 
ratiative  parameter,  is  prescrited  in  the  MABL  model  as  0.1. 
Long  wave  radiative  flux  is  calculated  differently  for 
the   cloudy   and    cloud    free   cases.  In    the    cloudy   case,      the 

only  clouds  permitted  are  non-black  stratus.  Radiative  flux 
is  a  function  of  cloud  emissivity,  e^,,  which  depends  on  the 
total  liquid  water  content,  W,  cf  the  cloud.  W  profiles  are 
nearly  linear  with  height  above  the  LCL  (Davidson,  et  al. , 
1984)    and,    along    with    e_  can   be   described    by: 


c 


W    =    0.5p^(h   -    Z^)g,^  (2.9) 

€^,=    1    -    exp  (-aW)  (2.  10) 

where  p  is  the  density  of  air  (1.25  X  10"3  gm/cm^)  ,  h  is 
the  height  of  the  mixed  layer  (cloud  top),  Z^.  is  the  LCL 
(cloud    bottom),    a    =    0.158    (Slingo,    et   al.,     1982),    and    g,^     is 

I  h 

the    liquid    water    content   at    the   cloud    top. 

Using  the  Stef an-Eoltzmann  law,  the  net  long  wave  cloud 
top  radiation  flux,  1 ph'  ^^^  ^^  calculated  from  the  cloud 
top  temperature,  T  ^  ,  and  the  effective  radiative  sky 
temperature,    T^.        ,    cloud    bottom   net   radiative   flux,    L„    ,    is 

'•  'sky'  'nc' 

calculated    in    the      same    way    using    the      cloud   bottom   tempera- 
ture,   T      ,    and   the    SSI,    T    .       These    fluxes    are   given    by: 

Lnh=  e^a(T    *    -    T    .   *)  (2.  11) 

nhc'h  sky'  *' 

1-00=  €^a(T3*    -    T^*)  (2.12) 
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where  a  is  Stefan's  constant  (4.61  X  10"ii)  and  €^  is 
obtained  from  e«iuation  2.10  .  Ihe  net  long  wave  flux  at  the 
surface,    F,       ,    becomes: 

long 

p         =    CT  [T    *    -  €^  T*   -     (1    -  €^)  I,.      *  ]  (2.  13) 

long  i-      s  c  *  c'      sky       -•  *  ' 

where  T  is  the  average  temperature  of  the  cloud- 

For  the  cloud  free  case,  the  net  fluxes  are  calculated 
at  the  top  of  the  mixed  layer  {Z  =  h)  and  at  the  surface  (Z 
=  0)  by  integrating  the  flux  emissivity  profile  (Fleagle  and 
Businger,  1980) ,  which  is  a  function  of  the  water  vapor  and 
temperature  profiles.  The  net  long  wave  flux  at  the  surface 
for  the  clear  sky  case  is: 

F,    =  F  -  F .  (2,  14) 

long      u     d  ^  ' 

where  F^  and  F^  are  the  upward  and  downward  radiative  fluxes 
respectively. 

C.       OCEANIC    BOONDABY    LAIEB     (OBI)    MODEL 

The  OBL  model  is  a  one-dimensional,  second  order  turbu- 
lence closure,  vertically  integrated  (bulk)  model  of  the 
upper  ocean  surface  turbulent  boundary  layer  or  mixed  layer 
developed  by  Garwood  (1977) .  It  uses  the  continuity  equa- 
tion for  an  incompressible  fluid,  the  first  law  of  thermody- 
namics (heat  equation),  the  conservation  of  salt  equation, 
an  analytical  equation  of  state,  the  Navier-Stokes  equation 
of  motion  with  the  geostrophic  component  eliminated,  and  a 
two-component   vertically    integrated    TKE    budget. 

An  understanding  of  the  dynamics  of  the  entrainment 
process  is  necessary  to  predict  the  MLD  change  with  time. 
The  stable  water  mass  underneath  the  mixed  layer  is  destabi- 
lized  and   eroded    by    TKE   in   the    aixed    layer.      This    TKE    budget 
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is  the  basis  for  an  entrainment  hypothesis.  A  closed  system 
of  equations  is  obtained  by  using  the  bulk  buoyancy  and 
momentum  ec[uations  with  the  mean  turbulent  field  modeling  of 
the  vertically  integrated  equations  for  the  individual  TKE 
components. 

To   better  define   the      mixing  process,      separate   vertical 
and    horizontal    equations      for   TFE    are    used.  Buoyancy   flux 

and  shear  production  provide  energy  for  vertical  mixing, 
with  buoyancy  flux  being  a  somewhat  more  efficient  source  of 
energy  for  mixing  due  to  its  direct  contribution  to  the 
vertical  component  of  TKE.  The  buoyancy  equation  is  derived 
from  the  heat  and  salt  equations  using  a  linearized  equation 
of   state: 

p    =  pj  1   -  «  ( 9   -    9^  }    -   /3  (S    -    S^)  ]  (2.  15) 

Buoyancy  is  given  by: 

b  =  g(p^-  '^)/p^  (2.16) 

where: 

9  =  temperature 

S  =  salinity 

g  =  gravity 

P  =  density 

«  =  thermal  expansion  coefficient 

/3  =  density  coefficient  for  salt 

The  tilde  represents  the  total  instantaneous  value  and  the 
subscript  "o"  denotes  an  arbitrary,  but  representative, 
constant  value.     In  the  short   term,   salinity   has  little 
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effect  on  the  upper  ocean  density  profile  except  at  higher 
latitiudes  (Miller,  1976)  .  Tliis  leaves  temperature  as  the 
dominant  density  factor.  However,  by  using  buoyancy  (b) 
instead  of  only  temperature  (9)  allows  the  model  to  be 
applied  in  cases  where  evaporation  and/or  precipitation 
contribute  significantly  to  the  surface  buoyancy  flux. 
Model  initialization  requires: 

(1)  mixed  layer  temperature  and  salinity  profiles; 

(2)  wind-driven  horizontal  current  profiles. 

The  initial  salinity  profile  is  made  isohaline  if  unknown 
and  the  initial  currents  are  assumed  to  be  zero  if  unknown, 
lack  of  initial  information  about  current  and  salinity 
profiles  is  not  a  serious  deficiency  because  the  model  will 
evolve  reasonable  transient  profiles  of  salinity  and 
momentum  after  only  a  single  diurnal  cycle,  and  the  final 
results  are  not  sensitive  tc  the  initial  current  and 
salinity  profiles  (Davidson  and  Garwood,  1984).  It  is  also 
possible  to  prescribe  an  upper  ocean  internal  vertical 
velocity  (upwelling  and  downwelling)  if  it  is  known  to  be 
significant.  The  mixed  layer  depth,  h,  is  defined  as  the 
shallowest  depth  at  which  the  observed  density  value,  a^,  is 
0.02  o^  units  greater  than  the  observed  surface  density 
value. 

At  each  one  hour  model  time  step,  the  following  boundary 
conditions  are  required: 

•  net   upward  turbulent   heat  flux   at   the  water   surface 
(sensible  and  latent  heat  flux)  plus  back  radiation; 

•  incident  solar  radiation; 

•  the  fraction  of  short  wave  radiation  absorbed  in  the  top 
one  meter  of  the  ccean; 

•  surface  wind  speed  and  direction; 
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••    clcud  cover ;  . 

•  temperature     (SST/        dry    bult   and    dew-point      air    tempera- 
ture) ;    and 

•  precipitation    (P). 

Vertical  velocity  at  the  bottoir  of  the  mixed  layer  may  be 
prescribed  hourly.  Other  physical  and  model  constants  that 
need  to  be  prescribed  are: 

•  the  extinction   coefficient  for   short   wave   radiation 
absorption; 

•  the  Coriolis  parameter; 

•  the  critical  Richardson  number   for  stability  adjustment 
below  the  mixed  layer; 

•  expansion  coefficient  for  temperature;  and 

•  density  coefficient  for  salt. 

The  model  forecast  is  not  particularly  sensitive   to  these 

constants,   but   reasonable  values   have  been   determined  by 

Gallacher,   et  al.   (1933),  and  they  need  not  be  readjusted 

for  geographic  and  seasonal  variability. 

By   use  of   bulk   aerodynamic   formulas,   the   turbulent 

fluxes  of  latent  heat,  Qg  /  and  sensible  heat,  Q^  ,  can  be 
estimated  as  follows: 

2e=  C,  {.98E3  -  EJO^^  (2.17) 

Qh=  C,(T^  -  Tja^^  (2.18) 

Net  back  radiation  is  estimated   from  the  empirical  equation 
(Husby  and  Seckel,  1578): 
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Q^    =  1.14  X  10-7(273.16  +  TJ*{.39  -  .SE^i/Z)        (2.19) 
(1  -  .6C2) 

where: 

E  =  saturated   vapor   pressur€  of  the   marine   air   (0.98 
corrects  for  salt  defects) 

£3=  vapor  pressure  of  air  bas€d  on  dew  point  temperature 

Tg=  air  temperature 

Tg =  sea  surface  temperature 

C  =  fractional  cloud  cover 

C^ =  drag  coefficient 

The  upward  heat  flux,  Q  ,  is  then  given  as: 

and  the  solar  radiation,  Q^,    by: 

Qs  =  (1  -  aa''  )  (1  -  .66C3)Q^  (2.21) 

The  constants  "a"  and  "b"  are  adapted  from  Tabata  (1964)  and 
the  cubic  cloud  cover  correction  from  Laevastu  (1960). 
Clear  sky  radiation,  Q^,  is  given  by  Seckel  and  Beaudry 
(1973): 

Qq  =  A^j  ■«•  A^  cos  (j)  +  B-,  sin  (j)  ■*■    A^cos    2(j)  (2.22) 

+  B2sin  20 

where  the  coefficients  (A^,  A^  ,  etc.)  were  calculated  by 
harmonic  representation   of   the   values   predicted   in   the 
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Smithsonian  Meteorological   Tables    (List    1958)    with: 

(f)    =     (2V365)  (t    -    2  1)  (2.  23) 

where  "t"  is  the  Julian  day  of  the  year. 

The  ocean  is  not  very  transparent  to  solar  radiation. 
In  the  open  ocean,  approximately  50%  is  absorbed  within  the 
first  meter  (this  fraction  is  called  "Rf").  Absorption 
varies  from  region  to  region  and  is  highly  dependent  on  the 
concentration  of  absorbing  particles  such  as  phytoplankton, 
"Gelbstoff"  (yellow  substance)  ,  and  suspended  particulate 
matter.  Coastal  regions  absorb  more  radiation  because  of 
increased  amounts  of  suspended  particulates.  Very  little  of 
this  absorbed  radiation  penetrates  below  the  mixed  layer 
because  of  the  low  thermal  conductivity  of  the  underlying 
stable  thermocline.  Thus  most  cf  this  energy  is  transferred 
upward  out  of  the  ocean  and  back  into  the  atmosphere.  What 
short  wave  radiation  remains  penetrates  the  mixed  layer  and 
is  attenuated  in  an  exponential  fashion  depending  on  water 
turbidity.   Net  heat  flux  at  the  surface  is  given  as: 

Qnet  =  C^  +  (Rf)Qs  -  Qs  (2.24) 

Surface  buoyancy  (heat  and  salt)  and  momentum  fluxes  can 
be  computed  by  use  of  the  foregcing  equations.  Mixed  layer 
turbulent  temperature,  salinity,  velocity  and  buoyancy 
fluxes  are  given  by: 


(T'W)  =  -Q^^^/qC  (2.25) 


(S'wM  =  (P  -  E)S„  (2.26) 


(u«w»)  =  u*2  (2.27) 


(b'w')  =  g[a(T'w»)  -/3(S»w»)]  (2.23) 
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where  the  subscript  "o"  refers  to  the  surface  value.  The 
friction  velocity   in   air,    u* ,    is  calculated   by: 

^*    =    ("s/Pa)  *^  (2.30) 

where  r  is  the  surface  stress  (dynes/cmS) .  The  fluxes  of 
momentum^  radiation,  latent  and  sensible  heat  at  the  sea 
surface  determine  shallowing  (retreat)  or  deepening  of  the 
mixed  layer  by  entrainment.  If  there  is  a  positive  buoyancy 
flux  (Q„g^<0  and  E>?)  the  MLD  will  deepen.  This  will  be  the 
case  at  night  or  whenever  long  wave  radiative  cooling  plus 
the  upward  turbulent  fluxes  of  heat  and  moisture  exceed 
solar  radiation.  Negative  buojancy  flux  will  result  in  the 
shallowing  of  the  MIC  and  is  caused  by  the  domination  of 
daytime  solar  heating  (net  warming)  at  the  surface,  provided 
wind    stress  is   not   large. 

A  schematic  of  the  ocean  mcdel  computations  is  shown  in 
rig.  2.5.  New  ocean  temperature,  salinity  and  wind-driven 
current  profiles  are  predicted  at  one-hour  intervals  and  are 
used    to   predict    the    new    MLD. 

D-       COUPLED    BOOHDAET    LAYER    MODEI 

Coupling  of  the  atmospheric  and  oceanic  models  was  first 
accomplished      by    O'Loughlin       (1982)  by    matching      momentum, 

sensible  heat,  latent  heat,  and  radiation  at  the  air-sea 
interface.  The      atmospheric      part   of      the      coupled      model 

computes  the  radiative,  heat  and  momentum  fluxes  which  are 
inputs  to  the  ocean  model.  Since  the  ocean  model  has  an 
hour  time  step,  and  the  atmospheric  model  has  a  half-hour 
time  step,  the  ocean  model  was  inserted  as  a  subroutine  into 
the    atmospheric    model  and   is      called   every    other    atmospheric 
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time  step.  Before  the  atmospheric  fluxes  can  be  utilized  by 
the  ocean  model,  they  must  be  converted  to  the  appropriate 
units  for  seawater- 

The  ocean  model  prediction  requires  the  vector  surface 
stress  for  the  ocean  turbulent  velocity  flux,  ^*w^*  This 
requires  that  the  wind  field  be  broken  into  horizontal 
components: 
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O^Q^   =   -(Sin    9)Uio  (2.31) 

O^Q^  =   -(cos   9)0,0  ^2.32) 

where  U^^  is  the  wind  speed  at  10  aeters  and  U^q^  ,  U^Oy  ^^^ 
the  east-west  and  north-south  horizontal  components  respec- 
tively. The  wind  direction,  9,  is  relative  to  true  north. 
Then    u*^2   ±s  given   by: 

T3  =  P3U*2    =    p^u*^2  (2.33) 

the  density  of  seawater,  p  ,  to  be  about  1  gm/cm''.  Finally, 
the  horizontal  components  of  the  turbulent  velocity  fluxes, 
u*^jj  2  and  u*^  2^  are  computed  hy. 


U*2  =  C^O.O'  =  C^t^lOx^  *  OlOy  ^)  (2.34) 

^*wx^  =  C^^iox  ^  (2.35) 

^*wy"  =  C^"lOy"  (2.36) 

SST  is  passed  from  the  HABL  mcdei  back  to  the  atmospheric 
model  for  use  in  the  next  time  step.  Fig.  2.6  shows  the 
interrelationships  between  the  atmospheric  and  oceanic 
models  in  the  coupled  model. 

The  coupling  process  has  the  potential  of  significantly 
enhancing  the  OBL  prediction  qualities  because  of  the 
improvement  in  the  boundary  conditions  inherent  in  a  coupled 
OBL-MABL  system.  Longer  period  reasonable  forecasts  should 
be  obtained  due  to  the  feedback  between  the  adjacent  turbu- 
lent boundary  layers  and  the  associated  thermodynamic  and 
dynamic  adjustments  in  each  layer.  Formation  of  stratus  or 
fog  in  the  MABL  model  is  particularly  important.   Accurately 
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INITIAL  PROFILES: 
ATM-  TEMPERATURE 
HUMIDITY 

OCN-  TEMPERATURE 
SALINITY 


Fig.  2.6    Flow  Chart  for  Ccupled  OBL  and  MABL  Model 

prescribed  incoming  solar  radiation  is  imperitave  for  a 
useable  CBL  forecast.  This  can  only  be  done  by  use  of  an 
MABL  model  such  as  has  been  described  in  section  B  (Davidson 
and  Garwood,  1 984) . 
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E.       D*      ESTIHATION       FBOM      TOfiBDIENT         KINETIC      ENEBGY        (TKE) 
DISSIPATION 

Direct  measurement  of  u*  on  a  ship  is  rather  difficult. 
One  popular  technique  which  is  used  extensively  at  the  Naval 
Postgraduate  School  is  the  dissipation  method.  An  excellent 
review  of  this  method  as  weU  as  the  bulk  method  is 
contained  in      Schacher,      et    al.         (1981).  The    dissipation 

method      involves   the      use    of      the      turbulent   kinetic      energy 
(TKE)       equation      which,      assumirg    steady   conditions      and    the 
divergence    term      to   be     small   relative      to    the      ether   terms, 
becomes: 

u*    =    {[€kz]/[  0^(2/1)    -    (z/L)  ]}0  (2.37) 

The  denominator  on  tie  right  hard  side  of  equation  2.37  has 
teen  parameterized,  where  the  ouantity  (l)^(z/L)  -  (z/L)  is 
substituted  by   (j)      following    Wyngaard   and   Cote    (1971): 

0^      =       (1    +    0.5    jz/Li^)3/2  for   z/L    <    0  (2.38) 

(l)g      =       (1    +    2.5     (z/1)  ^)  ^  for    z/L    >   0  (2.39) 

Determination  of  the  dissipation  rate,  6,  necessary  for  this 
technique  requires  that  the  inertial  subrange  follows  a 
Kolmogorov  spectrum,  i.e.,  where: 

S^j  (k*)  =  a  €  2/3 k* "5/3  .  (2,40) 

S^j  (k*)  in  equation  2.40  is  the  spectral  density  of  the  hori- 
zontal windspeed,  «  is  a  coefficient  with  a  magnitude  of 
0.52  (Champagne,  et  al- ,  1977),  and  k*  is  the  wavenumber. 
Equation  2-40  has  been  experimeitally  verified.  On  the  R/V 
Acania    a      hot    film   sensor   was      used    to    obtain   a      spectrum   of 


34 


high  frequency  velocity  fluctuations  from  which  the  TKE 
dissipation  u*  was  calculated-  It  is  important  to  note  here 
that  poor  velocity  turbulence  spectra  are  obtained  during 
periods  of  fog,  rain  cr  drizzle.  This  leads  to  the  dissipa- 
tion method  being  invalid.  Errors  may  also  occur  if  the 
relative  wind  is  less  than  3  m/s  allowing  condensation  to 
form  on  the  hot  film. 

Applying  Taylor's  hypothesis^  i.e.,  k*  =  2^f/u',  equation 
2.U0  may  be  rearranged  to  become; 

€  =  240  [S^(f)  ]^  (f/IT)  ^2  (2.41) 

where  f  is  frequency  measured  in  Hz.    The  inertial  subrange 
spectral  densities  are  integrated  between  two  wavenumbers,  k^ 
and  k^  respectively,   in  equation   2.40  following  Khalsa  and 
Bussinger  (1977)  : 


^Au^    =      /     S     rk)dk    =     (3/2)  cc  e2/3(k/2/3    -    x-0)  (2.42) 


Applying  Taylor's  hypothesis,  equation  2.42  now  becomes: 

a^j2  =  0.230  (€0")  ^  (f/^  -  f^"^)  (2.43) 

Rearranging  equation  2.43  and  combining  with  equations  2.37, 
2.38    and   2-39,    with    a  von    Karman   constant    of   0.4,    we    obtain: 

u*    =    2.81    CT^,      (z/[U(l)^(z/L)  ]}  0  {f/^    -    f -^}  1/2  (2.44) 

Substituting  f  =  50Hz,  f  =  5Hz,  the  measurement  height,  z  = 
20.0m,  equation  2.44  reduces  to  its  operational  form  as 
follows: 


U*  =  3.69(7^   f  "20m^e^^/^^  ^"'^  K'^'^^) 
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Because  z/L  is  a  function  of  the  friction  velocity,  u*, 
according  to  equations  2.38  and  2.39,  the  solution  to  equa- 
tion 2.45  involves  an  iterative  process  (Large,  1979). 
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III.  SYNOPTIC  SITDATION 

A-   MIXED  LAYEH  DYNAIICS  EXPEEIHENT  (MILDEX) 

The  Mixed  Layer  Dynamics  Experiment  (MILDEX)  was 
conducted  during  the  period  of  24  October  to  10  November, 
1983  near  34°N  and  126°W.  MILDEX  was  a  multi-group  and 
multi-platform  experiment  designed  to: 

•  provide  magnitudes  of  air-sea  energy  exchange  rates  for 
use  as  boundary  conditions  in  mixed  layer  modelling. 

•  evaluate  the  drag  coefficieit  from  TKE  dissipation  meas- 
urements and  its  dependence  on  swell  amplitude  and 
direction. 

•  to  provide  a  time  series,  from  radiosonde  measurements, 
of  MABL  structure  for  model  verification. 

•  provide  an  intercomparison  between  radiative  transfer 
algorithms  developed  at  Scripps  Institute  of 
Oceanography  and  measured  data. 

•  evaluate  the  utility  of  SODAR  as  a  ship-borne 
instrument. 

•  evaluate  radiative  transfer  models  of  the  atmosphere 
with  cloudiness  as  a  unigue  parameter. 

Meteorology  data  were  collected  on  the  R/V  Acania  as  well  as 
two  other  platforms  during  MILDEX.  R/V  Acania  data  collec- 
tion consisted  of  measurements  of  windspeed,  temperature, 
humidity,  radiation  and  atmospheric  pressure  in  the  surface 
layer,  and  sea  surface  temperature  (measured  by  thermistor), 
and  radiosonde  derived  profiles  of  atmospheric  windsfeeds 
and  directions,  tenperatures  and  hamidity.  Radiosonde 
launches  were  coordinated  with  satellite  pass  times. 
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Fig.  3,1  depicts  the  locations  on  the  bow  of  the  E/V 
Acania  of  mean  and  turbulent  leasurement  sensors.  There 
were  two  levels  of  these  measurement  sensors:  bow  mast  at  5 
meters  and  main  mast  at  20  meters.  In  addition,  there  were 
hourly  observations  of  sea  state  and  cloudiness. 
Conductivity,  temperature  and  depth  (CTD)  casts  were  also 
made  from  the  R/V  Acania.  Two  ether  vessels,  the  R/?  Wecoma 
and  R/P  FLIP  collected  meteorological  and  oceanic  data  near 
the  same  site  as  R/V  Acania.  During  the  periods  of  interest 
in  this  thesis,  the  R/P  FLIP  was  always  within  two  kilo- 
meters of  the  R/Y  Acania.  The  meteorological  data  has  been 
compiled  in  the  MILDEX  report  (Geernaert,  et  al. ,  1983)  . 
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Fig.  3.1    Schematic  of  the  R/V  Acania 
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Tiine  series  of  wind  speeds,  directions,  SST,  air  temper- 
ature, and  surface  pressure  in  Fig.  3.2  show  seme  of  the 
variability  experienced  during  KILDEX.  Frontal  passages  are 
evident  on  29  October,  31  Octcber,  and  06  November.  For 
most  of  the  experiment,  the  SST  was  greater  than  the  air 
temperature  (except  en  29  October  and  09  November) ,  so  the 
surface  layer  was  unstable  or  near  neutral-  Only  for  6-8 
hour  periods  on  29  October  and  09  November  was  the  surface 
layer  stable. 

Three  periods  (cases)  were  considered  for  purposes  of 
this  study.  These  cases  were  chosen  on  the  basis  of  strong 
wind  forcing  of  the  OBL  and  the  availability  of  reliable 
radiosonde  data  for  model  initialization-  Fig,  3.3  shows 
the  movement  of  the  R/V  Acania  for  all  three  cases.  The  S/V 
Acania  was  located  near  34, ON,  125. 4W  during  all  periods. 
All  times  are  Pacific  Daylight  Time  unless  otherwise  speci- 
fied, MABL  turbulent  kinetic  energy  dissipation  rate 
derived  u*»s  were  obtained  froE  velocity  fluctuation  meas- 
urements from  the  R/V  Acania  main  mast  (at  20  m)  and  a  bow 
mast  (at  5  m)  instruments.  MLE  data  were  obtained  from  CTD 
casts  from  the  R/V  Acania  for  Case  I  and  by  the  R/P  FLIP  for 
Cases  II  and  III. 

Conditions  associated  with  these  cases  was  dominated  by 
increasing  clouds,  swell  and  winds  with  intermittent  to 
heavy  rainfall.  Heavy  rain  showers  occurred  on  29  October 
when  the  first  of  several  fronts  passed  through  the  area. 
Since  the  presence  of  fog,  rain  or  drizzle  affects  the 
validity  of  the  hot  film  signals  used  in  the  dissipation 
technique,  care  was  taken  to  exclude  affected  results  during 
the  periods. 

Ocean  conditions  were  also  variable.  MILDEX  experiment 
participants  reported  large  scale  advection  and  upwelling 
occurred  during  the  course  of  the  experiment.  Significant 
internal  wave  activity  was  also  evident  and  is  being 
analyzed  by  other  investigators, 
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Fig.  3.2   Tine  Series  of  MILDEX  Rind  Speed  and  Direction, 
SST,  Air  Temperature,  and  Surface  Pressure.   Times  are  PDT. 
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Fig.    3-3        Movenent  of   R/Y    Acania:      Case   I    (solid) , 
Case   II    (dash)  ,    Case   III    (dot)  . 


B-       CASE   I     (0900,    31    OCTOBER   -    C900,    01    N07EMBEE) 

During  this  period,  the  B/7  Acania  experienced  weak  ccld 
frontal  passage  at  aiout  1700  oe  31  October.  Fig.  3.4  shows 
the  National  Meteorological  Center  (NMC)  surface  analysis 
for  Case   I.  The   synoptic    situation   was    that      a   low   center 
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was  located  in  the  northeast  Pacific  ocean,  and  a  high  doiii- 
nated  the  Pacific  north  of  Hawaii.  The  entire  period  had 
nearly  100%  cloud  cover  and  intermittent  rain  or  drizzle. 
Several  rain  squalls  were  noted  after  the  frontal  passage. 
Clouds  varied  from  altostratus  prior  to  the  frontal  passage 
to  stratus  for  the  rest  of  the  Case  I  period.  No  fog  was 
reported,  so  a  considerable  amount  of  hot  film  data  was 
availai)l€.   Fig.  3.5  is  a  copy  cf  GOES  WEST  visual  imagery- 

ks  shown  in  Fig.  3.2,  winds  were  initially  from  the  SSW 
with  a  shift  to  the  west  after  the  frontal  passage  with 
speeds  varying  from  7  to  9  m/s  ahead  of  the  front  and  from  3 
to  5  m/s  behind  the  front.  There  was  a  steady  decrease  in 
surface  pressure  from  1018  mb  to  1014  mb.  Swells  increased 
from  4  to  9  ft  after  the  frontal  passage  and  gradually 
decreased  to  5  ft  towards  the  end  of  the  observational 
period.  Air  temperature  dropped  by  1.5  C  after  the  frontal 
passage.  Throughout  the  period,  the  SST  was  greater  than 
the  air  temperature  with  the  greatest  difference  at  night  as 
expected.  Air-sea  temperature  differences  varied  from  -.2° 
to  -2.3  C.  The  MABL  surface  layer  was  near  neutral  to 
unstable  during  the  period. 

C.   CASE  II  (1600^  03  N07EHBEE  -  1600,  04  NOVEMBER) 

During  this  period,  a  stronger  cold  front  than  in  Case  I 
was  approaching-  Fig.  3,6  is  the  NWC  surface  analysis  for 
Case  II.  A  low  center  was  positioned  in  the  Gulf  of  Alaska, 
and  a  high  pressure  center  had  intensified  west  of  Baja 
California  and  south  of  the  Acania.  No  frontal  passage  was 
observed  for  this  case.  Winds  were  generally  from  the  NNW 
shifting  to  the  west  in  the  middle  of  the  observational 
period  and  then  returned  to  the  NNW  at  the  end  while  wind 
speeds  fluctuated  between  3  and  7  m/s.  The  surface  pressure 
was  very  steady  changing  only  by  2   mb  from  1021  to  1023  mb. 
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Swells  were  significaEt  and  froo  the  NW  and  ranged  from  7  to 
9  ft.  Fluctuations  in  air  temperature  never  exceeded  1-0  C. 
As  in  Case  I,  the  SSI  was  greater  than  the  air  temperature 
with  the  air-sea  temperature  difference  varying  from  -.6°to 
-1.6 °C.  In  this  case,  the  MAEL  surface  layer  was  unstable 
throughout  the  entire  period. 

Fig.  3.7  shows  GOES  WEST  visual  imagery  for  this  case. 
Cloud  cover  was  heavy  in  the  late  afternoon  and  again  at 
mid-morning  on  the  next  day.  Curing  the  night,  the  sky  was 
cloud-free.  Fog,  cumulus  and  stratus  clouds  were  evident  in 
the  afternoon.  Stratus  with  some  cirrus  was  present  the 
next  morning.  No  precipitation  was  recorded  during  the  Case 
II  period. 

D.   CASE  III  (1600,  Oa  NOYEHBEB  -  1600,  05  NOVEMBEB) 

The  cold  front  approaching  during  Case  II  dissipated  at 
the  beginning  of  the  period  and  the  area  came  under  the 
influence  of  a  weak  high  pressure  south  of  the  R/V  Acania. 
Fig.  3.8  is  the  NMC  surface  analysis  for  Case  III.  Another 
low  pressure  system  is  to  the  NW  with  an  associated  cold 
front.  This  second  system  did  not  appear  to  influence  the 
weather  or  the  seas  in  this  case.  Winds  initially  were  from 
the  NNW  and  gradually  shifted  to  the  SSW  as  the  high  moved 
through  the  area.  Speeds  generally  decreased  from  about  7 
to  2  m/s,  and  air  temperature  fluctuated  no  more  than  .7  C. 
Surface  pressure  varied  by  2  mt  from  1022  to  1024  mb  while 
sea  swells  were  from  the  NW  at  8  to  10  ft.  As  in  the 
previous  two  cases,  the  3ST  was  greater  than  the  air  temper- 
ature with   the  air-sea  temperature  difference   varying  from 

O  O 

-1.1  to  -2.1  C.   Again,  the  MABI  was  unstable  for  the  entire 
period. 

This  last  case  had  the  least  cloud  cover  due  to  the 
influence  of  the  high  pressure  system-   Fig-   3- 9  shows  GOES 
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WEST    visual  imagery   for   Case   III.  At   the   beginniDg    of    the 

period,  50%  cloud  cover,  consisting  mainly  of  stratocuaulus, 
was    evident.  This   changed   to      cirrus   and    finally      a    clear 

night-time  sky.  On  the  morning  of  5  Nov.  there  was  30%  to 
50/S  cloud  cover  (mostly  cirrus  and  some  cumulus)  with 
clearing  in  the  afternoon.  Nc  precipitation  was  recorded 
for   this  period    which  was  good    for   hot  film   performance. 
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17.  DATA  AN I  BES0LT5 

A-   APPROACH  TO  FOBCIfiG  0*  INTEEPEETATION 

For  each  of  the  three  cases,  the  u*  associated  with  the 
model  predictions  was  compared  vith  the  TKE  dissipation  rate 
{€)  derived  u*  and  the  model  MLD  was  compared  with  the 
observed  MLD.  As  described  previously,  the  epsilon  (e) 
estimates  necessary  for  u*  calculations  by  this  method  were 
obtained  from  hot  film  measurements  on  the  bow  (5  m)  and 
main  masts  (20  m)  of  the  R/V  Acania  (Fig.  3.  1)  . 

The  MA3L  model  is  configured  for  ten  input  winds,  so  the 
ten  greatest  and  least  winds  were  chosen  from  each  24-hour 
period.  These  will  be  referred  to  in  the  discussion  as 
"peak"  winds.  Linear  interpolation  was  used  for  data  points 
between  these  ten  points.  However,  a  forecaster  at  sea 
would  not  predict  high  and  low  winds,  but  would  most  likely 
forecast  average  winds  over  the  period.  Therefore,  another 
model  run  utilizing  2.5  hour  wind  averaging  was  input  over 
the  same  periods  to  see  if  the  use  of  averaged  winds  would 
have  any  deleterious  effect  on  the  MLD  computations. 

The  purpose  of  these  comparisons  is  to  evaluate  atmos- 
pheric forcing  (u*)  in  relaticn  to  predicted  and  observed 
MLD's.  There  are  three  u*'s  used  in  this  study:  those 
estimated  by  bulk  formulations  using  peak  winds,  by  bulk 
formulations  using  averaged  winds,  and  by  dissipation  esti- 
mates. In  the  following  discussion,  the  first  two  sets  will 
be  referred  to  as  model  u*  (peak  or  averaged)  and  the  last 
set  as  dissipation  u*. 
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B.   CASE  I  (0900^  31  OCTOBEB  -  C900,  01  NOVEMBER) 

This  first  case  resulted  in  the  poorest  comparison  of 
predicted  and  observed  values  of  any  of  the  three  cases. 
Initialization  data  is  summari2ed  in  Table  .  I.  Table  II 
lists  the  2.5  hour  averaged  winds  that  were  used  for  subseq- 
uent runs.  In  Fig-  4.1,  a  comparison  of  the  two  wind  inputs 
is  shown.  Fig.  4.2  shows  a  comparison  of  u*'s.  A  diffi- 
culty in  comparing  model  and  instrument  derived  values  for 
this  case  was  instrument  performance  and  equipment  break- 
down.  Both  of  these  problems  were  encountered. 


TABLE    I 

Case  One   Initialization 

Atmospheric  Values 

Pressure:    10  19,0    mb                    Sea    T 

emp    (C) :    18. 64 

Pot.    T(C):     17.9                               Spec. 
Lift.    Cond.    Lvl:    351.74m         Mvec 

Hum:     10.47 

tion:    0. 

Ocean    Values 

Mixed   Layer    Depth:    34.5m        Temp. 

Jump    (C)  :    -.  16 

Below  Layer    Gradient:    -.0155 

Winds 

Time(PDT)                   Speed  (m/s) 

Direction 

0900                                 8.3 

181.0 

1120                                  9.5 

176.0 

1345                                  7.4 

205.0 

1745                                  9.1 

199.0 

2030                                  3.4 

218.0 

2130                                  6.4 

215.0 

0140                                  3.5 

259.0 

0345                                  6.1 

247.0 

0505                                  2.4 

230.0 

0720                                  5.9 

294.0 

Eain,  drizzle,  and  low  relative  wind  speeds  affected  the 
accuracy  of  both   the  bow  and  main   mast  dissipation  derived 
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TABLE 

II 

Case 

One  Averaged 

iinds 

Time  (PDT) 

S£eed (m 

Z§) 

Direction 

0900 

8.3 

181.0 

1130 

8.9 

178.0 

1400 

8.2 

203.0 

1630 

8.6 

196.0 

1900 

4.6 

263.0 

2130 

5.4 

221.0 

0000 

5.7 

253.0 

0230 

4.5 

253.0 

0500 

4.2 

231.0 

0730 

4.9 

295.0 
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Fig.  4.1    Comparison  of  Case  I 
Initial  and  Averaged  Winds. 
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Fig-    4-2        Case    I   Priction   Velocities 

u*  values.  Post-processing  remcval  of  the  bad  data  resulted 
in  significant  data  gaps.  For  the  bow  mast,  dissipation  u* 
values  are  missing  from  15  16  tc  2315  and  from  060U  onward. 
Only  12  useable  observed  data  points  were  available  (as 
compared  to  48  for  the  model)  fcr  the  period.  The  field  log 
indicates  a  lack  of  confidence  in  main  mast  spectra  by  the 
observers  onboard  due  to  spurious  voltage  readings  and 
connector  problems,  so  data  from  it  were  not  included  at  all 
in    this    comparison. 

Fig.      U.2    shows    the    peak    wind    and   averaged    wind    model    u* 
and   the    dissipation    u*    values.  As   expected,       the    model    u* 

values  are  a  function  of  the  differences  in  the  input  winds, 
Where  the  averaged  winds  are  higher,  the  u*'s  are  higher. 
The    dissipation      u*    values    obtained      from    the   bow      mast    have 


54 


the  same  trend  but  are  lower  than  both  model  u*'s  in  the 
pre-frontal  period  (before  1700) .  They  are  close  to  the 
"peaked"  wind  model  predictions  for  the  rest  of  the  time. 
All  of  the  u*  dissipation  values  are  lower  than  the  averaged 
wind  model  predictions. 

These  results  of  model  and  dissipation  u*  values  were 
not  completely  expected  for  th€  frontal  regime.  Previous 
comparisons  of  bulk  and  dissipation  u*'s  in  frontal  regions 
(Geernaert,  et  al. ,  1985)  indicated  that  the  dissipation 
u*»s  were  higher  than  bulk  values  prior  to  a  frontal  passage 
and  then  the  magnitudes  reversed  after  the  passage.  Here 
the  dissipation  u*»s  remained  lower  than  the  bulk  value 
throughout.  Generally,  the  study  is  revealing  that  small 
pre-frontal  wind  speeds  (less  than  6  m/s)  produce  greater  u* 
deviations  (turbulent  vs  bulk) .  Additionally,  if  the  swell 
and  wind  oppose  each  other,  greater  u*  deviations  were  also 
evident.  These  were  minimized  if  the  wind  and  swell  were  at 
right  angles  to  each  other.  The  reasons  for  this  are  not 
entirely  clear.  However,  it  has  been  postulated  (Geenaert, 
et  al.,  1985)  that  the  direction  and  magnitude  of  the  swell 
relative  to  the  wind  (primarily  the  interaction  of  wind 
stress  with  short  wave  elements  on  long  wave  crests)  affects 
the  roughness  length,  z^,  which  in  turn  affects  u*  and  C^ . 

No  attempt  will  be  made  with  these  results  to  examine 
how  the  swell  interacted  with  the  wind  and  how  this  could 
have  modified  the  u*  values.  Perhaps  this  case  did  not  have 
the  same  kind  of  swell-wind  interaction  as  in  previous 
studies.   This  is  an  area  for  further  investigation. 

Fig.  4.3  compares  the  model  and  observed  MLD's.  The 
model  prediction  was  based  on  bulk  u*  values  shown  in  Fig. 
4.2.  In  this  case,  model  results  (solid  line)  early  in  the 
period  do  not  agree  with  the  observed  MID  (dashed  line) . 
While  both  curves  nearly  coincide  at  the  start,  the  model 
curve  shallows   much  more   quickly  although   it  reaches   the 
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Fig.  4.3   Case  I  Mixed  Layer  Depth 

proper  depth.  Then  the  model  deepens  rapidly  and  stays 
significantly  deeper  for  the  rest  of  the  case.  An  MID  error 
of  up  to  10  meters  can  Le  tactically  significant  for  naval 
forces  in  predicting  sonar  performance. 

The  early  shallowing  may  not  constitute  a  model 
prediction  error.  A  CTD  cast  would  give  the  thermal  MLD 
while  the  model  predicts  the  turbulent  boundary  layer  depth. 
It  takes  a  finite  amount  of  tiae  for  the  temperature  field 
to  adjust  to  the  turbulent  field.  The  expression  for  the 
lag  is: 


t         ~        [    {AT)U*3]/Q^2 


(U.l) 


where : 

T    =    temperature 
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t  =  lag  time 

Q^ =  surface  heating 

Thus  based  on  equation  4.1,  the  greater  the  u*,   the  greater 
the  lag  and  the  greater  the  heating  (Q^)  ,  the  less  the  lag. 

Using  2.5  hour  averaged  winds  showed  no  improvement  in 
MLD  prediction  (Fig-  4.3  dotted  line)  and  in  fact  gave  a 
poorer  prediction.  "While  the  two  wind  inputs  (Fig.  4.1)  are 
quite  different,  the  output  Mil's  are  quite  close  together 
except  at  1300,  This  is  a  pre-frontal  period  of  strong 
winds  and,  from  1200  to  1700,  the  averaged  wind  speeds  are 
higher  than  the  initial  winds.  Since  MLD  is  a  function  of 
the  wind  speed  cubed,  the  stronger  winds  would  give  deeper 
mixing.  This  is  again  true  for  the  period  from  2300  to  the 
end  of  the  run  where  the  averaged  winds  are  generally  higher 
than  the  initial  winds.  Overall,  one  can  conclude  that  for 
this  case  both  the  "peaked"  and  averaged  winds  give  poor  MLD 
predictions  and  these  could  have  adverse  effects  for  naval 
operations.  However,  if  the  dissipation  u*  values  for  0900 
to  1700  (which  are  lower  than  the  bulk  values)  had  been  used 
by  the  model,  the  model  MLD  would  have  been  closer  to  the 
observed. 

C.   CASE  II  (1600,  03  NOVEMBER  -  1600,  04  NOVEMBER) 

The  coupled  model  MLD  prediction  agreed  well  with  the 
observed  MLD  changes  in  this  case.  The  model  predicted 
unstable  ABL  conditions  throughout  the  period  which  agreed 
well  with  observations.  Table  III  contains  the  initializa- 
tion data  for  this  case  and  the  2. 5  hour  averaged  winds  are 
listed  in  Table  IV.  Fig.  4.4  compares  the  two  input  winds 
for  this  case. 

This  was  a  period  with  unreliable  dissipation  u*  data. 
Fig.    4.5  is   a  plot  of  the   various  friction   velocities. 
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TABLE    III 
Case  Two   Initialization 

Atmospheric  Values 

Pressure:    1022.4    mh  Sea    Temp    (C):    18.78 

Pot.    1(C):     17.1  Spec.    Hum:    9.32 

lift.    Cond.    Lvl:    476.71m         Aavection:    0. 

Ocean   Values 

nixed  layer   Depth:    12.0m        Temp.    Jump    (C)  :    -.2 
Below  layer   Gradient:    -.0094 

Winds 


Time  (PDT) 

Speed (m/s) 

Direction 

1600 

4.0 

040.0 

1745 

5.6 

310.0 

2120 

8.6 

036.0 

2230 

4.7 

319.0 

2300 

6.4 

355.0 

0135 

2.4 

270.0 

0400 

7.0 

281-0 

0705 

3.0 

350.0 

1005 

7.2 

008.0 

1420 

4.6 

336.0 

While  fog  or  rain  were  not  a  factor  in  dissipation  measure- 
ments, equipment  failure  was  a  problem.  This  led  to  a 
reduction  of  dissipation  u*  data.  The  bow  mast  data  were 
not  available  from  1700  to  2300,  2330  to  0217,  and  0900  to 
1600  when  a  major  equipment  repair  was  accomplished.  Bow 
dissipation  u*  data  from  12  thirty  minute  periods  were 
available.  Main  mast  dissipation  measurements  were  not 
available  from  1600  to  0752  and  from  0822  to  1340.  It  was 
also  absent  after  1410  and  only  data  from  4  thirty  minute 
periods  were  available. 

At  the  beginning  of  the  period,  the  bow  mast  dissipation 
u**s  are  close  to  the  bulk  deri^/ed  (model)  values-  However, 
during  the  middle  part  of  the  period,  the  bow  dissipation 
values  are  consistently  lower  than  the  bulJc  (model)  although 
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TABLE 

17 

Case  Two   Averaged 

Winds 

lime  (PDT ) 

S_geed  (_m 

IS) 

Direction 

1600 

3.9 

040.0 

1830      • 

4,5 

345.0 

2100 

5.7 

349.0 

2330 

4.3 

308.0 

0200 

3.6 

274.0 

0U30 

6.  1 

282.0 

0700 

4.6 

307.0 

0930 

6.6 

357.0 

1200 

6.8 

348.0 

1430 

5,4 

341.0 
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Fig.  4,4   Comparison  of  Case  II 
Initial  and  Averaged  Rinds, 
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Fig.  4,5   Case  II  friction  Velocity 

they  tend  to  have  the  same  time  variations.    The  two  bc¥  a* 
dissipation  values  near  0800  ar€  close  to  the  model  values. 

It  is  apparent  from  Fig.  U.6  that  the  model  predicted 
the  MID  well  enough  for  Naval  tactical  applications.  Ihe 
observed  MLD  (dashed  line)  shovs  some  differences  from  0200 
to  1000  which  are  probably  due  to  internal  wave  activity, 
advection  or  as  a  result  of  the  ship  moving  to  a  different 
area.  Observed  MLD  data  was  not  available  for  1100  to  1300. 
This  contributes  to  the  abrupt  MLD  rise  at  1400.  The  model 
(solid  line)  predicted  the  MLD  guite  well  especially  for  the 
deepening  part  of  the  period.  Better  agreement  between 
observed  and  model  shallowing  would  probably  be  evident  if 
MLD  data  would  have  been  available  for  1100  to  1300. 
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Fig-  4„6   Case  II  Mixed  Layer  Depth 

Comparison  of  the  peaked  winds  with  the  2.5  hour  aver- 
aged wind  MLD's  yields  some  interesting  results.  During  the 
deepening  phase,  the  averaged  winds  consistently  give  a 
shallower  MLD  (dotted  line)  than  the  initial  winds.  At  this 
time  the  "peaked"  winds  are  consistently  stronger  (Fig.  4.4) 
than  the  averaged  winds  (except  at  0100)  and  would  produce 
more  vigorous  mixing.  However,  during  the  shallowing 
portion,  the  averaged  wind  prediction  is  significantly 
deeper  than  the  peaked  wind  prediction.  This  time  the  aver- 
aged winds  are  predominantly  stronger  than  the  peaked  winds. 
Again  the  predicted  MLD  results  are  expected  since  the  MLD 
is  proportional  to  u*3.  If  the  dissipation  u*  values  for 
0200  -  0600  would  have  been  used  by  the  model,  the  predicted 
MLD  would  have  been  shallower  fcr  that  period  and  closer  to 
the  observed. 
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D.       CASE    III     (1600,     04    NO?EMBEE    -    1600,    05    N07EMBEB) 

As      in      the      previous      case,           the      coupled      model  MLD 

predictions   agreed    well   with      the    observed    conditions.  ABL 
conditions  were    unstaisle   throughout    the   case.      Table   V    lists 

the    initialization      data    utilized   for      this    model      run,  and 

Table      VI    lists      the    2.5      hour    averaged      winds.         Fig.  4.7 
compares   the   two  input  winds. 


TABLE   V 
Case  Three   Initialization 

Atmospheric  Values 

Pressure:    1022.3    mb  Sea.   Temp    (C):    18.70 

Pot-    T(C):     17.7  Spec.    Hum:    9.48 

Lift-    Cond.    Lvl:    508.0m  Advection:    0. 

Ocean    Values 

Mixed   Layer    Depth:    11.5m         Temp.    Jump    (C) :    -.1 
Below   Layer    Gradient:    -.0053 

Winds 

Speed  (m^) 

5.2 

6,5 
5.6 
4.5 
3.7 
5.2 
2.6 
2.  1 
4.8 
2.4 


Time  (PDT) 

1600 

1755 

2100 

2300 

0046 

0155 

0445 

0715 

0915 

1045 

Direction 

340, 

0 

338. 

0 

351. 

0 

352, 

0 

338. 

0 

325. 

0 

020. 

0 

296. 

0 

333. 

0 

271. 

0 

This  case  had  no  fog  related  problems  affecting  the 
dissipation  data  but  for  some  periods,  the  low  relative  wind 
speeds  degraded  dissipation  reliability.  Several  large 
dissipation  data  gaps  resulted.  For  the  bow  mast,  data  gaps 
exist  for  1600  to  00  18,  0100  to  0815,  0845  to  1200,  and  from 
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TABLE  YI 

Case 

Three  Averaged 

Hinds 

TimefPDT) 

S£eed  (m/s) 

Direction 

1830 

6.C 

340.0 

2100 

5.6 

351.0 

2330 

4.3 

346.0 

0200 

5.0 

3  54.0 

0430 

2.7 

004.0 

0700 

2.3 

311.0 

0930 

3.2 

272.0 

1200 

2.5 

264.0 

1U30 

2.3 

235.0 
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Fig.  4.7   Comparison  of  Case  III 
Initial  and  Averaged  iinds. 
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Fig.  4.8    Case  III  Friction  Velocities 

1331  onward.  Only  seven  reliable  dissipation  values  were 
calculated  for  the  bou  mast.  Data  was  not  available  for  the 
main  mast  from  1600  to  0018,  0  100  to  0315,  and  from  1410  to 
the  end  of  the  period.  This  resulted  in  four  thirty-minute 
dissipation  averages.  Fig-  4.3  compares  the  different  fric- 
tion velocities. 

For  the  entire  period,  the  peak  and  averaged  wind  model 
u*'s  are  very  close  due  to  the  small  wind  differences 
between  the  two  (Fig-  4.7).  Both  the  main  mast  and  bow  mast 
dissipation  u* ' s  agree  reasonably  with  a*'s  from  the  peaked 
wind  model  run.  This  is  particularly  true  in  the  second  half 
of  the  period.  Some  variation  is  expected  since  the  dissi- 
pation u*  comes  from  instantaneous  measurements,  while  the 
model  uses  interpolated  values  between  the  input  winds. 
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Fig.  U.9   Case  III  Mixed  Layer  Depth 

The  model  also  predicted  the  HLD's  well  (Fig.  U-9) 
although  there  is  greater  deviation  from  the  observed 
(dashed  line)  than  in  Case  II.  As  before,  the  actual  MLD 
shows  some  internal  wave  or  advective  features  (particularly 
around  0000  and  0400).  Good  observed  MLD  data  was  available 
for  the  entire  period. 

The  initial  model  run  (solid  line)  agrees  well  in  the 
initial  stages  but  deviates  during  the  above  mentioned 
internal  wave  or  advective  events.  The  sharp  "dip"  at  about 
0930,  not  present  in  the  model  cutput,  is  an  artifact  intro- 
duced by  use  of  a  spline  interpolation  in  the  plotting 
routine.  The  model  also  shallops  to  the  proper  depth  but  it 
occurs  about  one  hour  early. 
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Ccmparing  the  averaged  winds  MLD  with  the  peak  winds  MID 
shows  that  both  are  very  similar.  The  plot  of  input  winds 
(Fig.  4.7)  shows  that  both  the  initial  and  averaged  winds 
were  very  close  throughout  the  period  except  between  0800 
and  1000,  so  both  model  runs  should  produce  similar  MLD's. 
After  1100,  the  averaged  wind  MLD  (Fig.  4,9  dotted  line) 
shows  a  double  maxima  during  shallowing.  It  is  unknown  at 
this  time  why  the  model  predicted  such  maxima.  They  have 
been  slightly  exaggerated  by  use  of  the  spline  interpolation 
but  are  present  in  the  model  output.  As  with  the  initial 
wind  MID,  an  anomolous  "dip"  is  observed  at  1030  also  due  to 
the  spline  interpolation.  Both  the  initial  and  averaged 
wind  MlD's  are  within  operatioral  tolerances  for  Navy  use 
but  the  forecaster  should  be  well  versed  in  the  model  prop- 
erties and  limitations. 
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V.     CONCIDSIONS    AND    RECOMMENDATIOHS 

A.       CCMCIDSIOHS 

The  coupled  model  has  been  subject  to  various  past  anal- 
yses and  performed  well.  This  study  has  examined  the 
performance  of  the  model  in  regimes  that  stretch  its  capa- 
bilities.     Several   important   conclusions   can   be    drawn: 

(1)  The  model  overpredicts  clouds  which  affects  short  and 
long  wave  radiative  effects  on  boundary  layer 
dynamics. 

(2)  Good  comparisons  of  bulk  and  dissipation  u*  for 
unstable   conditions   are    made. 

(3)  In  near  neutral  or  slightly  unstable  regimes,  a* 
predictions  vary  systematically  probably  due  to  swell 
and   wind    stress   interaction. 

(4)  Better  predictions  are  achieved  if  peak  winds  instead 
of  averaged  winds  are  used.  However,  a  forecaster 
would  only  have  average  wind  prediction  estimates 
available  at  the  time  of  operationally  initializing 
the    model. 

(5)  Given  reasonable  fluxes  cf  short  and  long  wave  radia- 
tion, the  model  predicts  routinely  observed  varia- 
tions  in    ocean  mixed    layer    depth. 

(6)  It  is  a  potentially  powerful  tool  and  an  aid  in 
single-station  forecasting  to  a  Naval  Oceanographer 
who   is   aware    of   the    model   limitations. 

The  model  performed  very  well  within  its  limitations. 
However,  further  work  needs  to  be  done  before  it  can  be 
implemented   in    an    operational   rcle- 
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B.   EECCMMENDATIONS 

In  order  to  facilitate  operational  usage,  several  areas 
should  be  explored.  Some  of  the  aforementioned  problems 
could  be  alleviated  by  attempting  the  following 
recommendations: 

(1)  Continue  research  and  development  of  the  coupled 
model  with  direction  toward  operational 
iiplementation. 

(2)  When  clouds  form,  the  model  assumes  them  to  extent 
from  horizon  tc  horizon.  This  is  unrealistic  and  a 
statistical  probability  of  percentage  cloud  cover 
should  be  added.  Partial  cloud  cover  would  also 
alter  short  and  long  wave  effects  on  boundary  layer 
dynamics. 

(3)  Incorporate  swell  and  wind  interaction  effects  into 
u*  calculations  based  on  bulk  formulae. 

(4)  Continue  to  analyze  MILEEX  data,  particularly  for 
periods  where  there  are  more  complete  dissipation 
u*»s. 

(5)  Plan  for  model  coupling  to  other  predictive  models 
such  as  lEEPS  and  ICAPS. 

(6)  Involve  students  in  at-sea  data  gathering  experi- 
ments. This  would  enhance  understanding  of  thesis 
research  and  minimize  "false  starts"  due  to 
attempting  cases  with  faulty  or  limited  data. 

The  model  has  tremendous  potential.  Performing  further 
analysis  and  tweaking  the  model  will  result  in  a  powerful 
tool  for  the  fleet;  one  that  can  be  used  in  times  of  crisis 
when  extensive  satellite  networks  may  not  be  available  for 
use  by  forecasters. 
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